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Exact evidence for the spontaneous antiferromagnetic long-range order in the 
two-dimensional hybrid model of localized Ising spins and itinerant electrons 
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The generalized decoration-iteration transformation is adopted to treat exactly a hybrid model 
of doubly decorated two-dimensional lattices, which have localized Ising spins at their nodal lattice 
sites and itinerant electrons delocalized over pairs of decorating sites. Under the assumption of a 
half filling of each couple of the decorating sites, the investigated model system exhibits a remarkable 
spontaneous antiferromagnetic long-range order with an obvious quantum reduction of the staggered 
magnetization. It is shown that the critical temperature of the spontaneously long-range ordered 
quantum antiferromagnet displays an outstanding non-monotonic dependence on a ratio between 
the kinetic term and the Ising-type exchange interaction. 
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INTRODUCTION 



Exactly solvable lattice-statistical models are of appre- 
ciable scientific interest as they bring a valuable insight 
into diverse aspects of quantum, cooperative, and criti- 
cal phenomenafii^ii^i^i^i^iii^ It should be mentioned, how- 
ever, that sophisticated mathematical methods must be 
usually employed when searching for an exact treatment 
of even relatively simple interacting many-body systems, 
while an exact treatment of more realistic or more com- 
plex models is often quite unfeasible or is accompanied 
with a substantial increase of computational difhculties. 
From this point of view, the exact mapping technique 
based on generalized algebraic transformations^' be- 
longs to the simplest mathematical methods, which al- 
lows to obtain the exact solution of more complicated 
model from a precise mapping relationship with a simpler 
exactly solved model. Following Fisher's ideaS)^ an arbi- 
trary statistical-mechanical system (even of quantum na- 
ture) that merely interacts with either two or three outer 
Ising spins may be in principle replaced by the effective 
interactions between the outer Ising spins through appro- 
priately chosen decoration-iterationi^ or star-triangle^^ 
mapping transformations. Even although the concept 
based on the generalized algebraic transformations has 
been worked out by Fisher more than a half century agO;^ 
the algebraic mapping transformations were initially 
widely used to treat only lattice-statistical models con- 
sisting of the Ising spins (see Ref. [l3 and references cited 
therein) before this conceptually simple approach was fi- 
nally applied to hybrid models composed of the Ising and 
classical Heisenberg spinS j^^'^^i^^'^^i^^'^^ as well as, the 



Another interesting application of the algebraic 
mapping transformations has recently been suggested 
by Pereira et al^'^ when applying the generalized 
decoration-iteration transformation to an intriguing 
diamond-chain model of interacting spin-electron system. 
In this diamond-chain model, the nodal lattice sites are 
occupied by the localized Ising spins and the mobile elec- 
trons can freely move on a couple of interstitial decorat- 
ing sites symmetrically placed in between two localized 
Ising spins. The main aim of this work is to treat exactly 
an analogous two-dimensional (2D) hybrid model defined 
on doubly decorated planar lattices, which should provide 
a deeper insight into how itinerant character of the mo- 
bile electrons will influence phase transitions and critical 
phenomena of this interacting spin-electron system. 

The outline of this paper is as follows. In Section [III 
we first provide a rather detailed description of the model 
under investigation together with the most crucial steps 
of the exact mapping method, which enables us to ob- 
tain exact closed-form expressions for the critical tem- 
perature, order parameter and other relevant thermody- 
namic quantities. The most interesting results for the 
ground state, the finite-temperature phase diagram and 
thermodynamics are presented and detailed discussed in 
Section [ml Finally, the summary of the most important 
scientific achievement is mentioned with several conclud- 
ing remarks in Section [IVI 



II. 



MODEL AND METHOD 



Ising and quantum Heisenberg spins 
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Let us consider a hybrid lattice-statistical model of in- 
teracting spin-electron system on doubly decorated 2D 
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FIG. 1: A part of the doubly decorated square lattice, which 
has one localized Ising spin at each nodal lattice site (full 
circle) and two mobile electrons at each couple of decorating 
sites (empty circles). The ellipse demarcates the localized 
Ising spins and itinerant electrons described through the kih 
bond Hamiltonian 



Here, 4^_^ and c^^ ,^ (a = 1,2, 7 =t,i) denote usual 
creation and annihilation fermionic operators and a^^ 
is the standard spin- 1/2 operator with the eigenvalues 
(T^Q, = ±1/2. The hopping parameter t takes into account 
kinetic energy of the mobile electrons and the exchange 
integral J describes the Ising-type interaction between 
the itinerant electrons and their nearest Ising neighbors. 

A crucial step of our approach lies in the calculation 
of the partition function. Owing to a validity of the 
commutation relation between different bond Hamilto- 
nians [Ti.i,'Hj] = 0, the total partition function Z can be 
partially factorized into a product of the bond partition 
functions Zk 

Nq/2 Nq/2 

^ = E n Ti-fcexp (-/37ife) = E n (2) 

{cTi} k=l {ai} k=l 



lattices, which have one localized Ising spin at each nodal 
lattice site and two delocalized mobile electrons at each 
couple of decorating sites. The magnetic structure of 
the model under investigation is schematically depicted 
in figure [T] on the particular example of the doubly dec- 
orated square lattice. For further convenience, the total 
Hamiltonian can be defined as a sum over bond Hamilto- 
nians H = J2k '^k, where each bond Hamiltonian Hk in- 
volves all the interaction terms of two itinerant electrons 
from the fcth bond of the doubly decorated 2D lattice 

T^k =-t (^c|,i |Cj,2 I -I- 4i,i'^fe2,i + cL^r'^fei^T + 42,i'^fci,i) 

^'2^ki (4-1, t"^*:!,! ~ 4i,i'^fcia) 



where /3 = l/{kBT), fee is Boltzmann's constant, T 
is the absolute temperature, N is the total number of 
the Ising spins (i.e. the nodal lattice sites) and q is 
their coordination number (i.e. the number of nearest- 
neighbor decorating sites). Next, the symbol J2{cr } 
notes a summation over all possible spin configurations 
of the localized Ising spins and the symbol Tr^ stands 
for a trace over degrees of freedom of two mobile elec- 
trons from the fcth couple of decorating sites. An ex- 
plicit form of the bond partition function Zk can be 
subsequently acquired by a direct diagonalization of the 
bond Hamiltonian ([ij. The matrix representation of 
the bond Hamiltonian Hk in the orthonormal basis of 
states IV'j) = {4i,t42,t|0)' 4i,i42,il0)' cIijcLal^)' 
cIijcL^tIO)' 4i,T4ialO)' 42,T42,ilO)} l^'t>cls the 
empty state) reads 



hki — hk2 
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(3) 



where we have defined two new parameters hki = J'^kil'^ 
and hk2 = J^kil'^ that include the Ising-type interaction 
between the mobile electrons and their nearest-neighbor 
localized Ising spins. It is noteworthy that the parame- 
ters hki and hk2 can alternatively be viewed also as gener- 
ally non-uniform effective field, which is produced by the 
localized Ising spins on the itinerant electrons situated 
at the nearest-neighbor decorating sites. The straight- 
forward diagonalization of the Hamiltonian matrix ([3]) 



yields six eigenenergies 

Ekl,k2=^{hkl + hk2)i Ek3,k4 — 

Ek^.k6^±\j {hki - hk2f + U\ (4) 

which can be further used for the relevant calculation of 
the bond partition function Zk- After tracing out the 
degrees of freedom of the itinerant electrons, the bond 
partition function Zk merely depends on spin states of 
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two localized Ising spins and, besides, its explicit form 
immediately implies a possibility of performing the gen- 
eralized decoration-iteration transformation^jiSiii 



and the honeycomb lattici 



,32 



Zk=Y^ exp {-(iEk^) = 2 + 2 cosh 

i=l 

2 



2 cosh 



Aexp(/3i?(7|i(7^2)- 
(5) 



The physical meaning of the mapping transformation ([5]) 
is to replace the bond partition function Zk by the equiv- 
alent expression, which would contain the effective pair 
interaction R between the localized Ising spins only. The 
mapping parameters A and R are unambiguously given 
by the 'self-consistency' condition of the algebraic trans- 
formation ISI, which must hold for any combination of 
spin states of two Ising spins a^^ and cr^g involved therein. 
It can be readily proved that the decoration-iteration 
transformation ([5]) indeed represents a set of two inde- 
pendent equations, which directly determine the map- 
ping parameters A and R 



A^iViV2) 



1/2 



pR^2ln{Vi/V2) 



(6) 



that are for the sake of brevity expressed in terms of the 
newly defined functions Vi and V2 



Vi=2 + 2cosh(/3J/2) + 2 cosh(2/3t). 



t/2=4 + 2cosh f3^/j^ + (4t)2/2 



(7) 



Substituting the algebraic transformation (O into the 
factorized form of the partition function which phys- 
ically corresponds to performing the decoration-iteration 
mapping at each bond of the doubly decorated 2D lat- 
tice, consequently leads to a simple mapping relation- 
ship between the partition function Z of the interacting 
spin-electron system on the doubly decorated 2D lattice 
and, respectively, the partition function Zim of the sim- 
ple spin- 1/2 Ising model on the corresponding undeco- 
rated lattice with the effective (temperature-dependent) 
nearest-neighbor interaction R 



(8) 



The mapping relation ((8|) essentially completes our exact 
calculation of the partition function Z, since the partition 
function of the nearest-neighbor spin- 1/2 Ising model has 
been precisely calculated for several 2D lattices (for re- 
views see Refs. T^laSfail). For brevity, let us therefore 
merely quote the respective results for the partition func- 
tions of the spin-1/2 Ising model on the square latticeiS, 



lim i-lnZiM = ln2 + -i^ / / ln[cosh^ (/3i?/2) 



JQ 

sinh {f3R/2) (cos + cos (f>)]d0d(l) 



-1 p27r p27T 

/ / ln[{l + cosh3(/3i?/2) 
"T^ Jo Jo 



(9) 



lim — In Zim = In 2 , 

N-*ooN 16 

- sinh^ iPR/2) [cos 61 -I- cos + cos(6l + </>)]}/2]d6ld(?!). (10) 

In what follows, the precise mapping relationship ([5]) be- 
tween the partition functions will be used as a starting 
point for performing a rather comprehensive analysis of 
the critical behavior, the order parameter, and basic ther- 
modynamic quantities. 



A. Critical condition 

In order to locate a critical point of the investigated 
spin-electron system, one may take advantage of the fact 
that the partition function always becomes non-analytic 
at a critical point. It can be easily understood from 
Eqs. ([l])-® that the mapping parameter A cannot cause 
a non-analytic behavior of the partition function Z and 
thus, the investigated spin-electron system becomes crit- 
ical if and only if the corresponding spin-1/2 Ising model 
with the effective coupling (3R becomes critical as well. 
Accordingly, the critical points of the investigated spin- 
electron system can be straightforwardly obtained from a 
comparison of the effective temperature-dependent cou- 
pling PR with the relevant critical point of the corre- 
sponding spin-1/2 Ising model on the undecorated lat- 
tice. However, the mathematical structure of the map- 
ping parameter R has another important consequences 
on a critical behavior. More specifically, one may easily 
prove from Eq. ([7]) a general validity of the inequality 
Vi < V2, which in compliance with the definition ([5]) 
implies the antiferromagnetic nature of the effective in- 
teraction (i? < 0) in the spin-1/2 Ising model on the cor- 
responding undecorated lattice. Owing to this fact, the 
interacting spin-electron system is always mapped on the 
spin-1/2 Ising model with the antiferromagnetic nearest- 
neighbor interaction, which has a non-zero critical tem- 
perature only on the 2D loose-packed lattices such as 
squarei? and honeycomb^^ii^^'^i^'^ lattices. As a mat- 
ter of fact, it is well known dictum that the spin-1/2 Ising 
model with the antiferromagnetic nearest-neighbor inter- 
action on close-packed lattices like triangula r^^'^^i'^^ and 
kagomei2i^ lattices does not exhibit a spontaneous long- 
range order at any finite temperature and consequently, 
it does not have any finite-temperature critical point that 
would correspond to the order-disorder transition. 

Bearing all this in mind, we will consider hereafter only 
the interacting spin-electron system on the doubly dec- 
orated 2D loose-packed lattices. It is worthwhile to re- 
mark that the critical temperature of the antiferromag- 
netic spin-1/2 Ising model on the loose-packed lattices is 
equal to the one of the ferromagnetic model and hence, 
the critical points for the spin-electron system on the 
doubly decorated square and honeycomb lattices readily 
follow from the conditions Pc\R\ = 2 ln(l -f \/2)^'^ and, re- 
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spectively, Pc\R\ = 21n(2 + V^)^ where Pc 
and Tc is the critical temperature. 



B. Staggered magnetization 



l/(fcBTc) Callen-Suzuki identityilii^iiai^i^ 



As a direct consequence of the mapping equivalence 
with the antiferromagnetic spin-1/2 Ising model on the 
relevant undecorated loose-packed lattice, one should 
anticipate predominantly antiferromagnetic character of 
the interacting spin-electron system on the doubly dec- 
orated 2D loose-packed lattice as well. Let us therefore 
calculate the staggered magnetization as the most com- 
mon order parameter inherent to the antiferromagnetic 
spin alignment. Using the exact mapping theorems devel- 
oped by Barry et aZ.j2^ii2iiilii^ the spontaneous staggered 
magnetization of the localized Ising spins can easily be 
calculated from the exact spin identity 



1 



1 



K2)b. = 



miM(/3i?)(ll) 



where the symbols (. . .)t^j and (. . denote standard 
canonical ensemble average performed within the inter- 
acting spin-electron model on the doubly decorated 2D 
lattice and, respectively, its equivalent spin-1/2 Ising 
model on the corresponding undecorated 2D lattice."'^ 
The exact spin identity pip furnishes an accurate proof 
that the staggered magnetization of the Ising sublattice 
in the interacting spin-electron model on the doubly deco- 
rated 2D lattice directly equals the staggered magnetiza- 
tion of the spin-1/2 Ising model on the corresponding un- 
decorated 2D lattice with the antiferromagnetic nearest- 
neighbor interaction R < 0. However, the spontaneous 
staggered magnetization of the antiferromagnetic spin- 
1/2 Ising model on the loose-packed 2D lattices precisely 
coincides with the spontaneous magnetization of the fer- 
romagnetic model, i.e. the quantity, which has exactly 
been determined for several 2D Ising models (see Ref. |4J 
and references cited therein). In this regard, the stag- 
gered magnetization of the antiferromagnetic spin-1/2 
Ising model on the square^ and honeycomb^ lattices 
for instance read 



1 



miM = 



16a;^ 



1 - 



(l-a;2)4_ 
(l-a;)3(l-a;2)3 



(square) 

1/8 



(12) 
(honeycomb) 



where x = exp{—/3R/2). The above formulas complete 
our calculation of the staggered magnetization of the lo- 
calized Ising spins when substituting the exact expression 
p2)) with the appropriately chosen effective coupling ([6])- 
d?]) into the exact spin identity (fTTI) . 

On the other hand, the staggered magnetization of the 
itinerant electrons delocalized over pairs of decorating 
sites can be calculated with the aid of the generalized 



t x\ _ / Trfe/(cj'„ ,^,Cj^„ .y)exp(-/3?T!fc)' 



't,j 



(13) 



where a — 1,2, 7 =tiii ^^'^ f is in principle an arbi- 
trary function of the creation and annihilation fermionic 
operators from the fcth bond Hamiltonian With the 
help of the exact identity (fT3)) . the spontaneous staggered 
magnetization of the itinerant electrons can be calculated 
from the expression 



1 



I 1 ( dZk 



fci 



dhki 



t.j 
(14) 



where Si 



^ka,i'^ka,i 



)/2 marks the zth 



component of the spin operator of the mobile electron. 
After a straightforward calculation based on the differ- 
ential operator technique exp(ad / dx + bd/dy)g{x,y) — 
g{x + a,y + 6)^1^ and the exact van der Waerden iden- 
tity exp(ccr^) = cosh(c/2) -I- 2a^ sinh(c/2), the staggered 
magnetization of the itinerant electrons can be related to 
the staggered magnetization of the localized Ising spins 
through the precise formula 



J 



sinh 



V^^^ + W2 + cosh fyj^TM 



(15) 



Since the exact expression for the staggered magnetiza- 
tion of the localized Ising spins is already known from 
Eqs. (|ll |) -(fT2 |) . the formula p5|) provides the relevant 
exact result for the staggered magnetization of the itin- 
erant electrons hopping between the decorating sites. 



C. Thermodynamics 

Before concluding this section, it is worthy of notice 
that several basic thermodynamic quantities can also be 
easily derived from the exact mapping equivalence ([8]) be- 
tween the partition functions Z and Zim- For instance, 
the Hclmholtz free energy F, the internal energy U, the 
entropy S, and the specific heat C, can directly be com- 
puted from the basic relations of thermodynamics and 
statistical physics such as 



F^-fceTlnZ, U 



d\nZ 
dp ' 



dF_ 5C/ 
' dT' ~ dT' 



III. 



RESULTS AND DISCUSSION 



Now, let us proceed to a discussion of the most inter- 
esting results obtained for the interacting spin-electron 
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system on the doubly decorated 2D lattices. Before do- 
ing this, however, the relations ([SI)-© might serve in 
evidence that the effective nearest-neighbor interaction 
R of the spin- 1/2 Ising model on the corresponding un- 
dccoratcd lattice is invariant under the transformation 
J —t — J. A change of the ferromagnetic Ising interac- 
tion J > to the antiferromagnetic one J < actually 
causes only a rather trivial change of the mutual spin ori- 
entation of the itinerant electrons and their nearest Ising 
neighbors. This observation would suggest that the crit- 
ical temperature as well as other thermodynamic quan- 
tities remain unchanged under the sign change J — > — J 
and thus, one may further consider the ferromagnetic in- 
teraction J > without loss of the generality. 



A. Ground state 

First, let us take a closer look at the ground-state be- 
havior. It is quite evident that the ground state will 
correspond to the lowest-energy eigenvalue of the bond 
Hamiltonian ([T]) , which can be obtained from the eigenen- 
ergies ^ by considering four available configurations of 
the Ising spins Uki and ak2 explicitly involved therein. It 
turns out that the lowest-energy eigenstate constitutes a 
peculiar four-sublattice quantum antiferromagnet, which 
can be characterized through the eigenvector 



Nq/2 

fc=l 




1 



1 



J 



2t 



(16) 



|0), 



where the product runs over all bonds of the doubly dec- 
orated 2D lattice, the former ket vector determines spin 
states of the localized Ising spins, and the latter one spin 
states of mobile electrons. Interestingly, one may also 
find a much simpler goniometric representation of the 
lowest-energy eigenstate |AF) by introducing the mix- 
ing angle through the definition tan2(/) = 4t/J, which 
yields for the particular case with J > 0— 



.t 



• 2i t t 



Nq/2 

|AF)= W IT,i>a,i,o-fc2 [C0S-(/)4i -^,.j;2a 
k=l 

+ sin</'cos0(4i ,fc]:.^^^ -H42,T42,i) l^)' (1'^) 

Altogether, the four-sublattice quantum antiferromag- 
net |AF) can be characterized by a perfect Neel order 
of the Ising spins situated at the nodal sites of some 
loose-packed 2D lattice, while the mobile electrons de- 
localized over its decorating sites rest in the entangled 
state composed of two intrinsic antiferromagnetic states 
c^-^ |C^2 j,|0)i '^Ii x"^l2 T 1^^' non-magnetic ionic 



majority microstatc W"^^k3v!l^^^^W 



kV r r 



minority microstates w"^^I)~C/^^w 



a,. 



^k\ 



FIG. 2: A schematic representation of the quantum 
entanglement of two antiferromagnetic and two non- 
magnetic ionic states, which occurs in |AF) on account 
of a virtual hopping process of the mobile electrons. 
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FIG. 3: The probability distribution for the four entan- 
gled microstates as a function of the ratio between the 
hopping term t and the Ising exchange constant J. 



states cj^j^ |cj^j^ Jo), |cj^2 x|0)- The quantum entan- 
glement of those four microstates arises from a virtual 
hopping process of the itinerant electrons, which is dia- 
grammatically illustrated in Fig. [2l The respective prob- 
ability distribution of the four entangled microstates is 
displayed in Fig. [3] as a function of a relative strength of 
the hopping term. 

Let us make a few comments on an origin of 
the remarkable four-sublattice quantum antiferromagnet 
emerging in |AF). The hopping process of both itiner- 
ant electrons energetically favors their antiparallel spin 
alignment and this antiferromagnetic correlation is sub- 
sequently mediated through the Ising-type exchange in- 
teraction J also on their two nearest-neighbor Ising spins 
(hence, the antiferromagnetic character of the effective 
interaction i? < between the Ising spins). If a relative 
strength of the hopping term is sufficiently weak, the 
Ising spins prefer the ferromagnetic alignment with re- 
spect to their nearest-neighbor mobile electrons as it can 
be clearly seen from the occurrence probability of the ma- 
jority microstate cj,j^ ^c^j xl*^) ^^^^ converges to px^x ^ 1 
in the limit tj J — s- 0. It is nevertheless worth mentioning 
that the occurrence probabilities of three minority mi- 
crostates monotonically increase upon strengthening the 
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tIJ 

FIG. 4: The dimensionless critical temperature as a function 
of the relative strength of the hopping term for the interact- 
ing spin-electron system on the doubly decorated square and 
honeycomb lattices. 

hopping term at the expense of the occurrence probabil- 
ity of the majority microstate until they asymptotically 
reach the same value — = p^i.o = Pxt.o ^1/4 
in the limit t/ J — > oo. 



B. Finite-temperature phase diagram 

Next, let us examine critical phenomena associated 
with finite-temperature phase transitions of the sponta- 
neously long-range ordered phase |AF). It is worthwhile 
to recall that the critical temperature can easily be ob- 
tained by solving numerically the critical conditions de- 
rived in the foregoing section. The finite-temperature 
phase diagram in the form of the critical temperature 
versus the kinetic term dependence is shown in Fig. [4] 
for the interacting spin-electron system on doubly deco- 
rated square and honeycomb lattices. It is quite obvious 
from this figure that the interacting spin-electron system 
on any loose-packed doubly decorated 2D lattice exhibits 
the same general trends in the relevant dependences of 
the critical temperature. The critical temperature ini- 
tially exhibits a relatively rapid increase from zero tem- 
perature with increasing the ratio between the hopping 
term t and the exchange constant J until it achieves its 
maximum value. The critical temperature then gradually 
decreases upon further increase of a relative strength of 
the kinetic energy before it again tends to zero tempera- 
ture in the other particular limit t/ J oo. In the limit 
t/ J — > 0, the observed zero critical temperature can be 
attributed to a localization of the itinerant electrons at 
particular decorating sites. In fact, there does not ex- 
ist any other interaction between the itinerant electrons 
within the tight-binding model described by the Hamilto- 
nian ^ except the effective interaction originating from 
their virtual hopping process. If the hopping process of 
the itinerant electrons is ignored, the interacting spin- 
electron system then effectively splits into independent 



fragments each of them having one central Ising spin cou- 
pled to the q localized electrons from its nearest-neighbor 
decorating sites. In the other particular limit t/J^ oo, 
the zero critical temperature results from the same prob- 
abilities of the four entangled microstates of |AF) as it 
has been already discussed by the ground-state analysis. 
Owing to this fact, the decorating sites occupied by the 
mobile electrons have non-magnetic character, which is 
compatible with a disappearance of the effective interac- 
tion R = between the localized Ising spins that occurs 
in the particular limit t/ J — > oo. 

C. Order parameter 

Temperature variations of both sublattice staggered 
magnetizations rrii and nie are depicted in Fig. [5] for 
three different values of a relative strength of the kinetic 
energy. The most obvious difference between the sub- 
lattice staggered magnetizations rrii and rrie pertinent to 
the Ising spins and itinerant electrons, respectively, lies in 
the quantum reduction of the latter staggered magnetiza- 
tion. The greater a relative strength of the hopping term 
is, the greater quantum reduction of the staggered mag- 
netization rrie can be observed in Fig. [5] in concordance 
with the relevant ground-state prediction of the lowest- 
energy eigenstate (jl6l) . as well as, the zero-temperature 
limit of the expression both yielding 

^ ' 2 + (4t)2 ^ ' 

On the other hand, the sublattice staggered magnetiza- 
tion rrii pertinent to the Ising spins always starts from 
its maximum possible value that implies a perfect Neel 
long-range order of the localized Ising spins. For com- 
pleteness, it is also worthy of notice that both sublattice 
staggered magnetizations tend to zero in a vicinity of the 
critical temperature with the critical exponent /3 = 1/8 
from the standard Ising universality class. 

D. Specific heat 

Finally, thermal dependences of the specific heat are 
plotted in Fig. [6] for the interacting spin-electron sys- 
tem on the doubly decorated square lattice at three dif- 
ferent relative strengths of the kinetic term t/J = 0.1, 
0.25, and 0.5. It is quite evident that the investigated 
model system exhibits the familiar logarithmic singular- 
ity at critical temperature of the order-disorder transi- 
tion, which provides another independent confirmation 
of the critical behavior from the standard Ising univer- 
sality class. In addition to the logarithmic divergence ob- 
served in a close vicinity of the critical point, the specific 
heat also displays a round Schottky-type maximum in 
the high-temperature tail of the specific heat curve. This 
maximum is well separated from the logarithmic singu- 
larity at relatively weak hopping terms (see solid line for 
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FIG. 5: Thermal dependences of both sublattice staggered 
magnetizations for the interacting spin-electron system on the 
doubly decorated square lattice at three different values of a 
relative strength of the hopping term t/J = 0.2, 0.4, and 1.0. 
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FIG. 6: Temperature variations of the specific heat for the in- 
teracting spin-electron system on the doubly decorated square 
lattice at three different values of a relative strength of the 
hopping term t/J — 0.1, 0.25, and 0.5. 



t/ J — 0.1), while it becomes superimposed on a logarith- 
mic divergence at moderate strengths of the kinetic term 
(see dotted line for t/J = 0.25). Last but not least, the 
round high-temperature maximum again separates from 
the logarithmic singularity upon further increase of a rel- 
ative strength of the hopping term. The round maximum 
then shifts towards higher temperatures and it becomes 
the broader, the stronger a relative strength of the kinetic 
term is (see dashed line for t/J — 0.5). 



sites. It has been shown that the ground state of the 
investigated model system defined on any doubly dec- 
orated loose-packed 2D lattice represents an interesting 
four-sublattice quantum antiferromagnetic phase. It is 
worth noticing, moreover, that this four-sublattice quan- 
tum antiferromagnet exhibits a remarkable combination 
of the spontaneous long-range order manifested through 
a non-trivial criticality at non-zero temperatures with 
obvious macroscopic features of quantum origin such 
as the quantum reduction of the staggered magnetiza- 
tion pertinent to the itinerant electrons. To the best 
of our knowledge, the model under investigation thus 
represents a rather rare example of the exactly solved 
model with such an intriguing combination of otherwise 
hardly compatible properties. To compare with, the in- 
teracting spin-electron system at a quarter filling (i.e. 
with one mobile electron per each couple of decorating 
sites) merely exhibits a classical ferromagnetic or fer- 
rimagnetic spontaneous long-range order depending on 
whether the ferromagnetic or antiferromagnetic interac- 
tion between the localized Ising spins and itinerant elec- 
trons is assumed?^ This would indicate that an existence 
of the four-sublattice quantum antiferromagnetic phase 
is closely related to a collective motion of electrons even 
if there does not exist any direct interaction between mo- 
bile electrons within our tight-binding model. 

Besides the purely academic interest in searching ex- 
actly solvable quantum spin models, the considered 
model system has also been suggested in order to bring 
insight into a magnetism of hybrid systems consisting of 
both localized spins as well as mobile electrons. Among 
the most famous magnetic materials, which obey this 
specific requirement, one could mention the magnetic 
metal SrCogO 11—1^1^1^ or the series of polymeric co- 
ordination compounds [Ru2(OOCiBu)4]3[M(CN)6] (M = 
Pq60j61^^ Cr^iS3). ^jjg ij^ti-gj. family of magnetic ma- 
terials, one and three unpaired electrons of the trivalent 
metal ions such as Fe^+ {S = 1/2) and Ct^+ {S = 3/2) 
are localized at the corners of a simple square lattice, 
whereas three unpaired electrons are delocalized over the 
mixed- valent dimeric unit Ru2^ residing each bond of the 
square latticci^ In this respect, the magnetic structure of 
this series of coordination compounds closely resembles 
the one of the suggested model system even if the elec- 
tronic structure of the mixed-valent dimeric unit Ruj^ 
would surely require more complex Hamiltonian in or- 
der to describe the double-exchange mechanism in the 
mixed-valent Ru2^ dimeric unit. In this direction will 
continue our further work. 



IV. CONCLUDING REMARKS 
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